NASA Technical Memorandum 81789 


NASA-TM-81789 19800011797 


WING/STORE FLUTTER WITH NONLINEAR PYLON STIFFNESS 


Robert N. Desmarais and Wilmer H. Reed^ III 


APRIL 1980 




7y 


ttol 







t\V^ 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 


LArJoLt:'; center 

LISRhRY, NASA 
HAMPTON. .VIRGINIA 



WING/STORE FLUTTER WITH NONLINEAR PYLON STIFFNESS 


★ ★★ 
Robert N. Desmarais and Wilmer H. Reed III 
NASA Langley Research Center 
Hampton, Virginia 23665 


Abstract 

Recent wind tunnel tests and analytical stud- 
ies show that a store mounted on a pylon with 
"soft" pitch stiffness provides substantial in- 
crease in flutter speed of fighter aircraft and 
reduces dependency of flutter on mass and inertia 
of the store. This concept, termed the decoupler 
pylon, utilizes a low-frequency control system to 
maintain pitch alignment of the store during maneu- 
vers and changing flight conditions. Under rapidly 
changing transient loads, however, the alignment 
control system may allow the store to momentarily 
bottom against a relatively stiff backup structure 
in which case the pylon stiffness acts as a hard- 
ening nonlinear spring. Such structural nonlinear- 
ities are known to affect not only the flutter 
speed but also the basic behavior of the 
instability. This paper examines the influence of 
pylon stiffness nonlinearities on the flutter 
characteristics of wing-mounted external stores. 

Introduction 

Analytical investigations of aeroelastic sys- 
tems are usually based on linear theory which as- 
sumes both the structural and aerodynamic proper- 
ties to be independent of the amplitudes of 
oscillation. Aircraft structures typically exhibit 
nonlinearities, however, such as backlash or kine- 
matic deflection limits in moving control surfaces 
and in the connecting structure between wing and 
external stores. Studies of flutter of wings with 
control surfaces that contain structural nonlinear- 
ities (references 1-5) have shown that nonlineari- 
ties can affect not only the flutter speed of the 
system but also the characteristics of flutter 
motion. Similar studies in reference 6 investigate 
the effects of control system nonlinearities, such 
as actuator force or deflection limits, on perform- 
ance of an active flutter suppression system. 
Whereas flutter of a linear system is characterized 
by an exponential growth of oscillation amplitude 
with time, flutter of a nonlinear system may be 
amplitude limited. On the other hand, a nonlinear 
system which is stable with respect to small dis- 
turbances may be unstable with respect to large 
ones. Interest in this particular problem stems 
from studies in reference 7 of a passive wing/store 
flutter suppression concept known as the "decoupler 
pylon." These studies and others have shown that 
a store mounted on a pylon with "soft" pitch stiff- 
ness can provide substantial increase in flutter 
speed and reduce the dependency of flutter on the 
mass and inertia of stores relative to that of 
"hard-mounted" stores. By decoupling the influence 
of store pitch inertia on wing torsion modes, the 
frequency separation between flutter-critical modes 
is increased and the flutter speed is increased 
as indicated in figure 1. 
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The decoupler pylon uses a low frequency con- 
trol system to compensate for changes in static 
deflections of the store that would otherwise exist 
during maneuvers or airsoeed changes. Depending on 
the time constant of the store alignment system and 
the rate of change of load, however, the store may 
deflect enough to exceed the linear range of the 
soft pitch spring and "bottom" against a relatively 
stiff back-up structure. Under such conditions, 
pylon pitch stiffness varies in a nonlinear manner 
which depends on both the static preload and 
oscillation amplitude of the store. 

Symbols 

K spring constant of linear soft soring 
equivalent linear spring constant of 
nonlinear spring 

M elastic restoring moment about store oitch 
axis 

M-| fundamental Fourier component of store pitch 
moment 

Mq static moment required to deflect store 

_ against hard spring 

M static preload moment 

N ratio of hard spring constant to soft soring 
constant 

V flutter speed 

^nom flutter speed of linear system with nominal - 
design, stiff pylon (soring constant = NK) 
time 

t time at which store contacts hard spring 
6 describing function (6 = K /K) 

* 

6 describing function when M/M^ = 1.0 

0 store pitch angle 

0, amplitude of sinusoidal store pitch 
oscillation 

0Q pitch angle at which store contacts hard 

_ spring 

0 static store pitch deflection due to preload 

(jo circular frequency 

Nonlinear System 

An idealized representation of the nonlinear 
pylon suspension system is shown in the upper part 
of figure 2. The store is suspended from the wing 
by a pivot located near the wing elastic axis. The 
store Ditch frequency is controlled by a soft linear 
spring. The soring stiffness K is chosen to make 
the uncoupled store pitch freouency somev/hat lower 
than the fundamental wing bending frequency when the 
Store is rigidlyjnounted (see reference 7). A 
static preload M is assumed to act on the store 
as a result of loads such as a high-g pitch-uo 
maneuver. When the preload exceeds a value Mq, the 
associated static pitch displacement causes the 
pylon to contact mechanical stops at 0 = Oq. The 
stiffness then increases by a factor N over that 
of the soft spring. This situation results in the 
load displacement curve shown in the lov/er part of 
figure 2. 
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In this naper the preload and resulting 
disDlacement are such that the discontinuity 
in the load displacement curve at negative 6 
values is never reached so only the positive 
discontinuity, Bq* is analyzed. The nonlinear 
effects of damping are ignored because studies 
in reference 7 showed that flutter is insensitive 
to pylon damping. 


Describing Function Analysis Technique 

The analysis technique used is the "describing 
function" or equivalent linearization method (see 
references 2 and 4). In oarticular, reference 4 
presents comparisons of the results of the 
describing function technique with results using 
a complete nonlinear solution of the same 
problem to illustrate the validity of the method. 


less than 0^, the equivalent spring constant is 

the same as that of the linear soft spring, 
i.e., 6 = 1.0. With increasing oscillation 
amplitude, the system begins to contact the hard 
spring when 0 = 0^ and the equivalent spring 

then stiffens as the amplitude increases. This 
transition between a linear soft spring and 
a nonlinear hardening spring occurs in figure 4 
at the oscillation amplitudes where the curves 
for constant preload break away from the 5 = 1.0 
line. Conversely, when preload exceeds M , the 

equivalent spring constant is the same as the 
linear hard spring (6 = 20) for small amplitude 
oscillations. However, as oscillation amplitudes 
increase and deflections enter the soft spring 
range, the stiffness is characterized by that of 
a nonlinear softening spring. 


The basis of the describing function method 
is to assume a sinusoidal displacement and then 
compute the load developed in the nonlinear 
soring. The soring load is then expanded into 
a Fourier series. The fundamental is retained, 
and higher harmonics, which are assumed to be 
negligible, are discarded. The spring constant 
of the equivalent linear spring is then determined 
by obtaining the ratio of the coefficient of the 
load fundamental to the displacement amplitude, 
i<e = M^/Op where 


0=0+0^ sin wt 


M = M + Z 
n=l 


sin nojt 
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bilinear spring is apparent in figure 4. When 
the preload exactly matches the "bottoming" load, 
M/Mq = 1.0, the equivalent spring constant is 

independent of oscillation amplitude. In other 
words, the system behaves as though the spring 
were linear for all oscillation amplitudes. The 
describing function corresponding to this transi- 
tion region between a hardening spring and a 
softening spring is designated' 6*, and its 
value depends on N, the ratio of the two spring 
constants. For the case shown in figure 4, 

N = 20 and 6* = 2.412. The variation of 6* 
with N is indicated in figure 5. 


Flutter Calculations 


In this type of analysis, for each preload 
and displacement amplitude the problem is linear, 
consequently a linear flutter analysis may be 
made. The nonlinearity only appears as d change 
in equivalent linear spring constant when the 
preload or displacement amplitude is changed. 

With reference to figure 3, the describing 
function Is computed as follows: The load M, 

developed in the nonlinear spring, is expressed 
in terms of the store oscillation amplitude 0^ 

and the time f at which the store contacts 
the hard spring. This load is integrated over 
one cycle to give an expression for the preload, 

FT (thezeroth Fourier coefficient). For each 
specified prelo*-:' FT, this expression is then 
solved for This information is then used 
to compute the first few Fourier coefficients. 

All the Fourier cosine coefficients integrate 
to zero because the load function is symmetric 
about the zeros of the cosine function. The 
first Fourier sine coefficient, , is used 

to compute the equivalent spring constant, 

Kg = M^/0p and the describing function, 6 = K^/K. 

The second Fourier sine coefficient, which has 
a very small magnitude, is used to assess the 
validity of the method. No other Fourier 
coefficients were computed. 

Figure 4 is a plot of describing function 
6 = K^/K versus amplitude for various preload 

moments. Note that when the preload is less than 
and the static plus dynamic deflection is 


The describing function method is incorporated 
into flutter calculations as follows. First, 
the flutter velocity is computed as a function 
of equivalent pylon pitch stiffness using any 
standard linear flutter analysis technique. For 
convenience, this flutter boundary is expressed 

in terms of nondimensional ratios V/V and 

nom 

6 = Kg/K where is the flutter speed of 

a linear system with nominal pylon pitch stiffness 
typical of current fighter design oractice (in 
this case, 20 times the soft pylon pitch stiffness). 
Then, the family of curves relating the describing 
function to amplitude and preload (figure 4), 
is cross plotted against the flutter boundary to 
eliminate 6. The result is a family of flutter 
boundaries (velocity versus displacement amplitude 
curves), one for each preload. This procedure 
is illustrated graphically in figure 6. The 
linear theory flutter boundary is the upper left 
part of the figure and the describing function 
(figure 4) is plotted below it. Flutter boundaries 
for the nonlinear system are in the upper right 
part and the dashed lines with arrows indicate 
how they are related. 

Several points can be made regarding the 
flutter boundaries shown in figure 6 for various 
preloads and oscillation amplitudes. Consider 
first the case where static deflections due to 
preload lies within the_l inear range of the 
decoupler pylon, i.e., M/M^ <1.0. If the system 

is initially at rest, the flutter onset speed 
is the same as that for the linear system, i.e.. 
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approximately V = 2.8 However, for speeds 

within the interval 2.3 < < 2.8, the 

system can experience divergent flutter oscilla- 
tions if a sufficiently large disturbance causes 
oscillations into the destabilizing stiff-spring 
range. For example, in figure 6 with a oreload 
of rVM = 0.6, it can be seen that the degrading 
effects^of stiffness nonlinearity on flutter 
begin to appear when the disturbance amplitude 
becomes greater than 0-|/9 q = 0*^ 

V > 2.8 Vnom, divergent oscillations occur for any 
disturbance. The line separating initial disturb- 
ances which cause flutter from those vihich do not 
is shown dashed. Note also that when M/Mq = 1.0, 
the flutter velocity is independent of the 
magnitude of the disturbance, the same as for a 
linear system. This interesting feature is, of 
course, a consequence of the descriMng function 
being independent of amplitude at M/Mq = 1.0 
for a bilinear spring, as was discussed earlier. 

Consider next, the range of preloads which 
cause static deflections greater than the limits 
of the soft spring, i.e., M/Mq > 1.0. Flutter 
onset under these conditions occurs at 

the flutter velocity of the linear system with 
stiff pylon. However, unlike a linear system, 
where oscillation amplitude grows indefinitely 
and exponentially with time, the flutter ampli- 
tude for this nonlinear system is self limiting 
because as amplitude increases the resultant 
softening effect of the nonlinear spring (see 
figure 4) tends to stabilize the system. This 
can be illustrated, for example, by the flutter 
point shown in figure 6 where, for 

the flutter amplitude is limited to about 0.1 6^. 

The solid curve shows the magnitude of the 
limited amolitude flutter. 

By physical reasoning it can be deduced that 
the sign of the slope of the velocity versus 
deflection amplitude curves, dV/d9], determines 
whether flutter oscillations will be divergent 
or of limited amplitude: a curve with positive 

slope indicates limited amplitude flutter; a 
curve with negative slope indicates the distur- 
bance amplitude that must be exceeded to cause 
divergent flutter. 


Applications 

In this section of the paper some analytically 
predicted effects of pylon stiffness nonlinear- 
ities on wing/sotre flutter are presented for 
two configurations: the F-16 and a flutter 

research model. Wind tunnel tests of both 
configurations with linear decoupler pylons have 
been conducted in the Langley Transonic Dynamics 
Tunnel. Additional tests with nonlinear pylon 
stiffness are being planned for the research 
model . 

F-16 Flutter Model 

The F-16 store configuration considered in 
this example is designated configuration 32 in 
reference 6. It consists of a GBU-8B store 
carried at wing stations 120, and an'AIM-9 
missile at each wing tip. The GBU-8B stores 


are mounted on decoupler pylons which give an 
uncoupled store pitch frequency of 4.0 Hz (on 
the full-scale airplane). This frequency, 
selected on the basis of a criterion of 
reference 7, is approximately 70% of the first 
antisymmetric wing bending frequency with the 
store rigidly attached. The pitch stiffness 
when the system bottoms against mechanical stops 
is taken to be the stiffness of the nominal 
F-16 pylon design which is 20 (N = 20) times 
greater than that assumed for the decoupler 
pylon. A linear flutter analysis for F-16 
configuration 32 was performed by General 
Dynamics, Fort Worth, with varying pylon pitch 
stiffnesses. Results of the analysis for 
anti symmetrical flutter (the most critical mode) 
at Mach number 0.9 are presented on the left 
side of figure 7. The companion curves plotted 
on the right side of figure 7 is a family of 
flutter boundaries which account for pylon 
stiffness nonlinearities. These flutter 
boundaries are for the same configuration as 
those previously presented, in abbreviated form, 
in figure 6 to illustrate the analysis procedure. 
Thus, earlier comments on figure 6 regarding 
divergent flutter and limited amplitude flutter 
are applicable to figure 7 as well. 

To put in better persoective the magnitude 
of static deflections a decoupler pylon with 
soft pitch stiffness might experience in pitch-uo 
maneuvers, some calculations have been made 
for the F-16 configuration considered in figure 7. 
In these calculations the benefits of an align- 
ment control system v/ere neglected and the 
most critical combinations of design maneuvers 
specified in Mil Spec MILA-8591E (ref. 8) v;ore 
assumed, i.e., pitch acceleration, +4 radians 
oer second; normal acceleration, +6 g's; and 
longitudinal offset of store center of gravity, 

3.5 inches forward. The static pitch deflection 
of the GBU-8B store predicted for this extreme 
pitch-up maneuver was only 0.9 degrees. 

Therefore, even though the decoupler pylon has 
a pitch stiffness that is low relative to the 
nominal pylon design stiffness, static deflection 
of the store during maneuvers is small. The 
bounds on store pitch deflection for other types 
of transient loads, such as gusts, are also 
being analyzed in a separate study. It appears, 
however, the main requirement for a store align- 
ment control system is to compensate for the 
quasi-steady variation in drag loads due to 
changing flight conditions. 

Decoupler Pylon Research Model 

The second configuration analyzed is the 
decoupler pylon research model used in studies 
reoorted in reference 7. The model is a 
cantilevered rectangular wing with a store 
mounted at the 80-percent semi scan station. The 
linear-system flutter (velocity) boundary for 
the model as a function of pylon pitch stiffness 
shown in figure 8 was derived from figure 6 of 
reference 7. As in the previous example, pylon 
stiffness is expressed in terms of the describing 
function for a bilinear spring with N = 20. 

Note that in contrast with the F-16 flutter 
boundary (figures 6 and 7) which decreases mono- 
tonically with increasing pylon pitch stiffness, 
the linear-system flutter boundary in figure 8 is 
characterized by a peak near 5 = 5,0. This peak 
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falls in the region of pylon stiffness where 
there is frequency coincidence between the 
uncoupled pylon pitch mode and the wing funda- 
mental bending mode. Because of this peak in ‘ 
the linear-system flutter boundary, the effects 
of stiffness nonlinearities on flutter are some- 
what different from, and more complicated than, 
the case shown oreviously. Again solid lines are 
used to show the magnitude of limited amplitude 
flutter and dashed lines lines to seoarate dis- 
turbance regions which lead to either stable 
motion, or catastrophic divergent flutter 
oscillations. For this nonlinear system, 
limited amplitude flutter can occur for any 
preload condition; iii the previous example, it 
occurred only v/hen M/M^ > 1.0. 

Conclusions 

This paper investigates the effects of pylon 
stiffness nonlinearities on the flutter charac- 
teristics of wings with externally mounted 
stores. In particular, the focus of the paper 
is on a passive wing/store flutter suppression 
concept known as the decoupler pylon. This i 

concept uses a soft pylon pitch spring to | 

decouple the store pitch mode from wing torsion 
modes assisted by a low frequency active control 
system to reduce static deflections of the store i 
due to maneuvers and changing flight conditions. j 
The structural nonlinearity under consideration 
is associated with bottoming of the system 
against a relative stiff backup structure as a 
result of excessive static and/or dynamic 
deflections of the store in pitch. By use of an j 
approximate analysis technique (describing ' 

function technique), the nonlinear flutter 
behavior of two wing configurations with external 
stores is studied. On the basis of these studies 
the following conclusions may be drawn: 

(1) If the store static pitch deflection due 
to preload falls within the linear stiffness 
range of the decoupler pylon (M/M^ < 1.0), the 

flutter speed is substantially greater than (more 
than twice) the flutter speed with a nominal 
stiff-pylon design. When store pitch oscilla- 
tions are superimposed on this static deflection, 
causing the system to bottom against a stiff 
backup structure, the flutter speed is changed 
but remains well above the flutter speed of the 
nominal stiff pylon. 

(2) If the stcre static pitch deflection 
due to preload exceeds the linear range of the 
soft pitch spring (M/M > 1.0), flutter onset 
occurs at the same speed as for the linear 
system with stiff pylon; however, the flutter 
amplitude is limited due to the stabilizing 
effect of the softening pitch spring with 
increasing deflection amplitude. 

(3) If the static preload exactly equals 
tjie load required to bottom the system 

(M/M^ =1), the flutter speed becomes independent 

of amplitude, as in a linear system, but is a 
function of N, the ratio of hard spring 
stiffness to soft spring stiffness. 


(4) Some sample calculations for an F-16 
fighter performing a design-limited pitch-up 
maneuver were made to determine the static pitch 
deflections of a decoupler-pylon mounted store. 
The maximum predicted deflection (without an 
alignment control system) was less than 1°. 

Thus, the requirement for an active alignment 
control system to avoid excessive static 
deflection of the store appears to be governed 
more by drag loads than by maneuver loads. 
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Fig. 1 Effect of decoupler pylon on 
flutter speed. 



Fig. 3 Representation of nonlinear spring 
by an equivalent linear spring. 





Fig. 4 Describing function for bilinear 
spring (N = 20). 


Fig. 2 Decoupler pylon with nonlinear 
pitch stiffness. 


5 


3.0 


M/Mq = 1.0 



Fig, 5 Describing function for transition 
between hardening spring and 
softening spring that occurs v;hen . 

M = M . 

0 



Fig. 7 Calculated flutter boundaries for 
1/4-scale F-16 model with nonlinear 
pylon stiffness (config. 32, ref. 7). 


Fig. 6 1 11 list ‘"at ion of method for determining 

wing/siore flutter boundaries with 
nonlinear pylon stiffness. N = 20. 




Fig. 8 Calculated flutter boundaries for 

flutter research model with nonlinear 
stiffness (config. 1, ref. 6). 
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